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CALCULATION  OF  THE  AERODYNAMIC  LOADING  OF  FLEXIBLE  WINGS 
OF  ARBITRARY  PLAN  FORM  AND  STIFFNESS 
By  Franklin  W.  Diederich 


SUMMARY 


A  method  is  presented  for  calculating  the  aerodynamic  loading,  the 
divergence  speed,  and  certain  stability  derivatives  of  wings  and  tail 
surfaces  of  arbitrary  plan  form  and  stiffness.  Provision  is  made  for 
using  either  stiffness  curves  and  root-rotation  constants  or  influence  , 
coefficients % in  the  analysis.  Computing  forms,  tables  of  numerical 
constants  required  in  the  analysis,  and  an  illustrative  example  are 
included  to  facilitate  calculations  by  means  of  the  method. 


INTRODUCTION 


The  distribution  of  the  aerodynamic  loading  on  wings  and  tail 
surfaces  is  important  both  for  the  structural  analysis  of  these  com¬ 
ponents,  since  it  determines  the  applied  bending  moment  and  torque 
acting  at  any  station,  and  for  their  aerodynamic  analysis,  since  it 
affects  the  stability  derivatives  to  a  large  extent .  At  high  speeds  the 
aerodynamic  loading,  particularly  in  the  case  of  swept  wings,  is 
greatly  affected  by  the  structural  deformations  caused  by  the  loading. 

The  present  paper  is  concerned  with  the  determination  of  the  effects  of 
structural  flexibility  on  the  aerodynamic  loading  of  wings  of  arbitrary 
plan  form  and  stiffness. 

The  problem  of  load  distribution  was  analyzed  for  unswept  flexible 
wings  as  early  as  192 6  (reference  l)  but  has  received  relatively  little 
attention  since  that  time.  The  only  new  effect  considered  in  subsequent 
work  is  aerodynamic  induction  (reference  2).  No  work  appears  to  have  been 
done  on  the  loading  of  flexible  swept  wings.  The  related  problem  of  aero- 
elastic  divergence  of  swept  wings  with  certain  prescribed  stiffness 
variations  has  been  treated  in  reference  3. 

I 

The  present  paper  treats  the  problem  of  aerodynamic  loading  by 
matrix  methods.  Aerodynamic  induction  is  taken  into  account  approxi¬ 
mately^  since  suitable  aerodynamic  influence  coefficients  are  not  avail¬ 
able  for  wings  of  arbitrary  plan  form.  When  they  become  available  they 
can  readily  be  incorporated  in  this  method.  Structural  flexibility  is 
taken  into  account  in  the  form  of  either  calculated  stiffness  variations 
or  measured  influence  coefficients.  The  required  integrating  matrices 
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are  presented  for  "both  a  six-point  and  a  ten-point  solution.  For  the 
six— point  solution  convenient  computing  forms  are  included  as  well.  The 
method  is  illustrated  by  means  of  an  example.  In  addition  to  the  analysis 
of  the  aerodynamic  loading,  the  determination  of  the  related  divergence 
speed  and  of  certain  stability  derivatives  is  discussed. 

For  the  convenience  of  the  reader  unfamiliar  with  matrix  terminology 
a  summary  of  matrix  methods  has  been  included  in  the  appendix.  The 
sections  on  "  Application  of  the  Method"  and,  in  particular,  "Instructions 
for  Solution"  may  be  rep,4  without  reference  to  the  section  "Derivation  of 
the  Method," 


SYMBOIS 


W 


a 

a* 

ac 

b 

c 

c 


ljw 


°BMw 


'W 


/me(lsA2el^cr2cos4A 


(GJ) 


aspect  ratio 
aeroelastic  matrix 
dimensionless  parameter  ^ 

parameter  (meqswcr  cos  a) 

section  aerodynamic  center,  measured  from  leading  edge, 
fraction  of  chord 

-wing  span,  inches 

chord  measured  parallel  to  the  air  stream,  inches 
average  -wing  chord,  inches 

/  2 1 

section  lift  coefficient 


qc 


/LwN 

wing  lift  coefficient  (  — 

VIs/ 


wing  root  bending-moment  coefficient 


.qSb/ 


wing  rolling-moment  coefficient  (2CSI^) 
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3 


[%] 

M1 

Wi 

[%] 

K] 

N 


m^qsy^c-r  tan  A  cos3a 


dimensionless  parameter 


"bending  stiffness  in  planes  perpendicular  to  the  elastic 
axis,  pound— inches2 

location  of  elastic  axis  measured  from  leading  edge,  fraction 
of  chord 

distance  from  reference  axis  to  section  aerodynamic  center 
(positive  forward ),  fraction  of  chord 

torsional  stiffness  in  planes  perpendicular  to  the  elastic 
axis,  pound— inches^ 

unit  matrix 

matrix  defined  by  equation  (12) 

integrating  matrix  for  single  integration  from  tip  to  root 
first  row  of  matrix 

integrating  matrix  for  double  integration  from  tip  to  root 
first  row  of  matrix 

integrating  matrix  for  single  integration  from  root  to  tip 

matrix  relating  concentrated  and  accumulated  torque 

matrix  relating  concentrated  loads  and  accumulated  bending 
moments 

matrix  converting  torques  due  to  distributed  loads  to  torques 
due  to  concentrated  torques 

matrix  converting  bending  moments  due  to  distributed  loads  to 
bending  moments  due  to  concentrated  loads 

lift  on  both  wings  but  excluding  lift  on  part  of  wing  covered 
by  fuselage,  pounds 

running  air  load  along  the  reference  axis,  pounds  per  inch 

accumulated  bending  moment  (in  planes  perpendicular  to  the 
reference  axis  unless  specified  otherwise),  inch— pounds 

effective  section  lift-curve  slope  for  angles  of  attack  due  to 
deformation,  per  radian 
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effective  section  lift-curve  slope  for  additional— type  angle— 
of-attack  distributions,  per  radian 

section  lift-curve  slope,  per  radian 

concentrated  load,  pounds 

wing-tip  helix  angle 

root— twist  constants  (see  equation  (9)) 
root-bending  constant  (see  equation  (9)) 
dynamic  pressure,  pounds  per  square  inch 
concentrated  torque,  inch-pounds 

total  wing  area  including  part  of  wing  covered  by  fuselage, 
square  inches 

distance  from  wing  root  to  wing  tip  perpendicular  to  the  air 
stream  (see  fig.  l),  inches 

length  of  wing  along  reference  axis  (see  fig.  l),  inches 

accumulated  torque  (in  planes  perpendicular  to  the  reference 
axis  unless  specified  otherwise),  inch-pounds 

distance  between  the  effective  root  and  the  innermost  complete 
section  of  the  torsion  box  perpendicular  to  the  elastic 
axis,  inches 

lateral  ordinate  measured  from  wing  root,  inches 
lateral  center  of  pressure,  incheB 
angle  of  attack,  radians 

equivalent  angle  of  attack,  radians  + 

local  dihedral  angle  due  to  deformation  or  slope  of  wing 
deflection  curve  at  reference  axis,  radians 

structural  deflection,  inches 

distance  along  reference  axis,  inches 

angle  of  sweepback  (measured  to  the  referenoe  axis  unless 
specified  otherwise),  degrees 
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N 

N 


influence-coefficient  matrix  for  wing  twist  in  planes  parallel 
to  the  air  stream  due  to  concentrated  unit  loads  applied  at 
the  reference  axis,  radians  per  pound 

influence-coefficient  matrix  for  wing  twist  in  planes  parallel 
to  the  air  stream  due  to  concentrated  unit  torques  applied 
in  planes  parallel  to  the  air  stream,  radians  per  inch- 
pound 

angle  of  twist  in  planes  perpendicular  to  the  reference  axis, 
radians 


Subscripts : 


a 

c 

D 

fw 

g 

LE 

M 

MAC 

P 

r 

rw 

s 

sub 

spr 

T 

TE 

w 


additional 

midchord 

divergence 

flexible  wing 

geometric 

leading  edge 

due  to  bending  moment 

pertaining  to  the  mean  aerodynamic  chord 

damping  in  roll 

at  root  or  effective  root 

rigid  wing 

structural  (due  to  structural  deformations) 

subsonic 

supersonic 

due  to  torque 

trailing  edge 

wing  exclusive  of  fuselage 
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[  ] 
[°] 


square  matrix 
diagonal  matrix 

DERIVATION  OF  THE  METHOD 
Method  Employing  Stiffness  Curves 


Assumptions.—  In  the  development  of  the  method  the  following 
'  assumptions  are  made: 

r 

(a)  The  effectB  of  aerodynamic  induction  may  he  taken  into  account 
by  applying  an  over— all  correction  to  strip  theory  and  rounding  off  the 
resulting  load  distribution  at  the  tip. 

(b)  All  deflections  and  angles  of  attack  are  small. 

(c)  The  wing  is  mounted  flexibly  at  an  effective  root  perpendicular 
to  the  elastic  axis  through  the  intersection  of  the  elastic  axis  and 
the  fuselage  (see  fig.  l),  the  root  rotations  being  proportional  to  the 
root  bending  moment  and  root  torque, 

(d)  An  elastic  axis  exists  in  the  outer  portion  of  the  wing,  this 
axis  being  defined  as  the  elastic  axis  the  wing  would  have  if  it  were 
mounted  rigidly  some  distance  outboard  of  the  root  approximately  perpen¬ 
dicular  to  the  midchord  line.  (Near  the  root  the  elastic  axis  is  defined 
as  the  extension  of  the  outboard  elastic  axis.) 

(e)  All  deformations  are  given  by  the  elementary  theories  of 
bending  and  of  torsion  about  the  reference  axis,  which  in  this  case  is 
the  elastic  axis. 

Air  loads.-  In  keeping  with  assumptions  (a)  and  (b)  the  force  on 
a  wing  section  of  unit  width  parallel  to  the  direction  of  flight  is 


2*  =  qc  ^rngOg  +  m^otg  j  cos  A 
=  mgqca  cos  A 


where  the  equivalent  angle  of  attack  a  is  defined  hy 


(1) 
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fl 
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•v 
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The  effective  section  lift-curve  slope  mg  for  angle— of -attack  distribu¬ 
tions  due  to  structural  deformations  has  been  given  in  reference  3  on 
the  basis  of  the  reasoning  of  reference  If-  for  subsonic  speeds  as 


“e  =  mo 


A 

A  +  4  cos  A 


(2) 


Similarly,  the  effective  lift— curve  slope  m^  for  additional— type 
angle— of -attack  distributions  is  determined  by  the  same  reasoning  as 


=  *0 


A 

A  +  2  cos  A 


(3) 


Both  slopes  must  be  multiplied  by  cos  A,  as  in  equation  (l),  in  order 
to  apply  to  loads  acting  on  sections  and  angles-  of  attack  measured 
in  planes  parallel  to  the  direction  of  flight. 


The  torque  of  the  running  load  V  about  the  reference  axis  is  l^pC 
for  uncambered  sections  (for  cambered  sections  the  torque  at  zero  lift 
must  be  added  and  the  analysis  of  the  following  paragraphs  modified 
accordingly).  .This  torque  may  be  resolved  into  a  running  torque  about 
the  elastic  axis  and  a  running  bending  moment  about  a  line  perpendicular 
to  the  elastic  axis.  The  running  load,  torque,  and  moment  must  then  be 
multiplied  by  cos  A  to  yield  their  values  per  unit  length  along  the 
elastic  axis,  so  that 

l  =  m^ca  cos2A  (4) 

or,  in  matrix  notation, 

}  =  meq(cos2A)  [c]|a  j  (4a) 


The  running  torque  and  the  running  bending  moment  are,  respectively, 

j^2©qC  cos  A]  and  j— Ze^c  sin  A  j>.  The  running  bending  moment  leads  to 

accumulated  bending  moments  which  have  to  be  added  to  the  accumulated 
bending  moment  due  to  the  running  load. 


The  accumulated  torque  T  is  obtained  from  the  running  torque  by 
an  integration  inboard  from  the  tip.  This  integration  may  be  performed 
by  a  matrix  |K-J  'which  is  based  on  Simpson*s  rule  with  a  modification 

suggested  by  T.  M.  Falkner  at  the  tip.  (See  appendix. )  The  effect  of 
Falkner*s  modification  is  to  round  off  the  calculated  load  distribution 
and  cause  it  to  go  to  zero  with  an  infinite  slope  at  the  tip,  as  the 
aerodynamic  lift  distributions  actually  do.  The  matrix  is  given  in 
table  I. 


8 


NACA  TN  No.  1876 


Similarly,  the  accumulated  bending  moment  M  is  obtained  by  a  double 
integration  inboard  from  the  tip  of  the  running  load  and  a  single 
integration  of  the  running  moment.  The  double  integration  may  be 
performed  by  another  matrix  [k2J  (given  in  table  II ),  which  iB  based 

on  the  equivalent  of  Simpson,s  rule  for  moments,  Faltaier's  modification 
again  being  made  at  the  tip.  The  derivation  of  the  integrating  matrices 
is  discussed  in  somewhat  greater  detail  in  the  appendix. 

The  accumulated  torque  and  bending  moment  may  then  be  written  as 


{t)  =  eA[%]  jle-jc  cos  a| 


=  me^sAelrcr2(GOs3A)[Kl]  (^)  {a} 


(5) 


and 


[m]  =  ba2[k^]  (2)  -  sa[k1]  {leqc  sin  a] 


me<lsA2cr 


N  (; 


-V 


®lr  \CrJ 


{a}  (6) 


*  Equations  of  equilibrium.-  The  equations  of  equilibrium  of  a 
deformed  wing  referred  to  the  elastic  axis  are 


GJ  §£  =  T  (7) 

dy 


El  —  =  M  (8) 

dy 

These  equations  must  be  integrated  outboard  from  the  root  to  obtain  <p 
and  T.  The  integrations  may  be  performed  by  a  matrix  [K3J  (see 

table  TTT  and  appendix),  also  based  on  Simpson* s  rule  without  the  tip 
modification,  however,  since  the  torques  and  moments  go  to  zero  with 
finite  and  zero  slopes,  respectively.  To  the  deformations  obtained  in 
this  manner  the  rotations  due  to  the  root  deflection,  cpr  and  Tr, 
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must  be  added.  The  root  rotations  are  defined  by  four  dimensionless 
constants: 


Sr  ~ 


W 


T. 


v/(GJ). 


(9a) 


Si 


v/(GJ). 


(9b) 


rrr/Tr 

Qr  «  — 

T  w/(El) 


(9c) 


<sr  -  -Ul 

M  w/(El) 


(9d) 


which  may  be  combined  into  two  other  constants 


Sr 


SA 


M=Tg5£ 


ctaVTr 

(s  ■ 

(gjL 

> 

\ 

w/(GJ)r 

'  (ex)j,  t8“  A%y 

J  cos  A 

(9e) 

%/^r  | 

/ 

^(>M  ' 

(gjL 

> 

v/(GJ)r 

-  7 - r—  tan  AQP%  J 

(Ei)r  rMy 

1  cos  A 

(9f) 

as  in  figure  1. 

The  deformations  may  then  be 

written  as 

h][(S,rj 

1  V 

+  — 

8a 

W  *  i  S,W{«}  I 

(10) 

If  equations  (5),  (6),  (10),  and.  (ll)  are  substituted  in  the  matrix 
equivalent  of  equation  (13),  the  following  relation  is  obtained.' 

{%}  =  a[A]{a|  ■  (l4) 

where  the  aeroelastic  matrix  [a]  is  defined  by 
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and 


d 

a 


(EE)r  ©os  A 


tan  A 


(l6a) 


are  substantially  the  same  parameters  as  those  used  in  reference  3  • 

Solution  of  the  equations.-  If  it  is  desired  to  calculate  the 
aerodynamic  loading  corresponding  to  a  given  geometrical  angle-of-attack 
distribution  and  dynamic  pressure,  equation  (l4 )  may  he  rewritten  as 
follows : 


►  m 

[&]  “  afA|^  5  >  =  {ag}  (17) 

In  this  form  it  constitutes  a  set  of  linear  simultaneous  equations  for 
the  a  values  in  terms  of  ag  values.  From  the  calculated  a  values 

the  lift  distribution  may  be  determined  from  either  equation  (l)  or  (4). 

The  divergence  dynamic  pressure  may  be  obtained  from  equation  (17) 
by  setting  the  determinant  of  the  square  matrix  on  the  left  side  of  the 
equation  equal  to  zero.  This  procedure  is  equivalent  to  setting  <xg 

equal  to  zero  in  the  term  a  of  equation  (l4),  so  that 

{<%}  =  a  [A]  {03}  (18) 

The  critical  value  of  a  is  then  determined  by  matrix  iteration  and 
hence  the  divergence  dynamic  pressure  from  equation  (l6). 


Method  Employing  Influence  Coefficients 

The  assumptions  made  in  the  preceding  sections  concerning  the 
behavior  of  the  wing  structure  are  unnecessary  if  influence  coefficients 
for  the  given  structure  are  available  from  test  data  or  refined  methods 
of  calculation.  The  coefficients  most  convenient  for  this  analysis  are 
those  giving  the  rotation  of  the  structure  in  planes  parallel  to  the 
direction  of  flight  due  to  vertical  loads  applied  along  a  convenient 
reference  axis  and  due  to  torques  about  lines  perpendicular  to  the 


12 


NACA  TN  No.  1876 


direction  of  flight..  Since  it  is  usually  more  convenient  to  apply 
concentrated  rather  than  distributed  loads  in  structural  tests,  the 
influence  coefficients  are  considered  in  this  analysis  to  have 
been  obtained  in  this  manner. 

The  angle  of  structural  deformation  a0  may  be  expressed  in  terms 
of  the  influence  coefficients  ^  and  $p  as  follows: 

{<%}  -  +  MM  (19). 

where  the  R’s  and  P’s  are  arbitrary  concentrated  torques  and  loads, 
the  latter  being  applied  at  the  reference  axis.  The  accumulated  torques 
and  bending  moments  about  lines  perpendicular  and  parallel,  respectively , 
to  the  direction  of  flight  may  be  related  to  the  concentrated  torques 
nnfl  loads  by  means  of  the  summation  matrices  and  pwjj  (see 

appendix)  as  follows: 

|t’|  =  [Ei|.]{E}  “  tan  A{M*}  (2°) 

{**’}  -  «»[K5]{p}  <21) 


These  relations  may  “be  solved  for  the  values  of  R  and  P  required  to 
produce  given  distributions  of  accumulated  torque  and  bending  moment 


(r}=  [kJ  1|{t»}  +  tan 

(22) 

1? 

M =  ^H_1W 

(23) 

The  accumulated  torques  and  bending  moments 
load  are  then 

produced  by  the  air 

{t'I  -  3„[kJ  /i.'ejtj  -  (m'J  tan  A 

(24) 

{m«}  - 

(25) 

NACA  TN  No.  1876 


13 


Upon  substituting  equations  (22),  (23),  (2b),  (25),  and  (l)  into 
equation  (19),  the  following  equation  is  obtained: 


|aB}=  a*[A*]|aj  (2  6) 

where 


a*  =  meqs^Cj.  cos  A  (27) 

and 

0  2“|  0 

M-  <28) 


where,  in  turn 


M  -  N_1[ki] 

(29a) 

_  -1 

N  -  py  psj 

(29b) 

are  given  in  tables  IV  and  J. 

The  solution  of  equation  (26)  is  obtained  in  the  manner  previously 
described  for  equation  (it). 


APPLICATION  OF  THE  METHOD 
Determination  of  the  Structural  Parameters 


At  the  time  an  aeroelastic  analysis  is  performed  no  exoerimental 
stiffness  data  are  usually  available,  so  that  the  calculated  stiffness 
curves  must  be  used.  In  order  to  use  these  curves  it  is  necessary  to 
assume  the  existence  of  a  reasonably  straight  elastic  axis.  The  location 
of  thiB  axis  may  be  estimated  by  considering  it  to  be  the  line  connecting 
the  shear  centers  of  the  individual  sections.  If  the  elastic  axis 
obtained  in  this  manner  is  not  reasonably  straight  within  a  few  percent 
of  the  chord,  the  results  of  the  analysis  may  not  be  sufficiently  reliable. 
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The  stiffnesses  GJ  and  El  do  not  have  much  physical  significance 
inboard  of  the  last  point  where  there  is  a  complete  cross  section  of  the 
torsion  hoi.  (See  fig.  1.)  In  order  to  arrive  at  estimates  of  the  root 
stiffnesses  (GJ)r  and  (Ej)r,  which  serve  primarily  as  reference  values 

in  this  analysis,  the  stiffness  curves  have  to  he  extended.  It  is  con¬ 
venient  to  consider  the  stiffnesses  to  be  constant  inboard  of  the  last 
complete  section  of  the  torsion  box;  this  procedure  should  yield  conserva¬ 
tive  values  of  the  root  rotations. 

The  most  difficult  problem  incurred  in  analyzing  the  deflections  on 
the  basis  of  stiffness  curves  appears  to  be  the  estimation  of  the  root 
rotations.  As  used  in  this  analysis,  they  are  the  torsion  and  bending 
deflections  imposed  by  the  triangular  inner  portion  of  the  wing  and  the 
carry— through  bay  on  the  rest  of  the  wing.  As  seen  in  figure  2,  which 
is  plotted  from  the  data  of  reference  5,  these  values  are  essentially 
constant  along  the  span,  so  that  they  actually  constitute  rigid— body 
rotations.  (The  bending  rotations  have  been  obtained  by  talcing  the 
difference  in  slope  between  curves  calculated  by  considering  the  wing 
to  be  cantilevered  at  the  effective  root  -  the  root  used  to  calculate 
torsional  deformations  in  reference  5  -  an!  "the  averages  of  the  leading- 
edge  and  trailing-edge  deflections  actually  measured.  The  twists  were 
obtained  by  subtracting  the  twists  calculated  on  the  basis  of  the  assumed 
effective  root  from  the  measured  twists. ) 

rTThiq  rotations  should  in  any  practical  case  be  calculated  by  analyzing 
the  triangular  root  and  the  carry— through  bay  and  made  dimensionless  by 
means  of  equations  (9 ) •  If  such  an  analysis  is  not  available,  the  dimension¬ 
less  rotation  parameters  shown  in  figure  2  may  be  used  as  a  guide;  it 
must  be  kept  in  mind,  however,  that  in  the  case  of  a  sweptforward  wing 
the  parameters  Qq^  and  would  have  the  opposite  sign. 

Once  the  structure  under  investigation  is  built,  fairly  simple 
deflection  tests,  similar  to  those  performed  in  reference  5,  may  be  used 
to  check  the  root-rotation  parameters  by  calculating  the  differences 
between  the  observed  rotations  and  those  calculated  by  simple  beam 
theory  considering  the  wing  cantilevered  at  the  effective  root;  at  the 
Rflinfl  time  the  existence  and  estimated  location  of  the  elastic  axis  may 
be  verified.  If  the  experimental  program  is  fairly  extensive  it  is 
desirable  to  measure  influence  coefficients  directly .  These  influence 
coefficients  can  then  be  used  in  conjunction  with  the  alternate  method 
described  in  the  preceding  section  to  obtain  a  quick  check  on  the  aero- 
elastic  analysis  based  on  calculated  stiffnesses. 

The  influence  coefficients  used  in  the  analysis  consist  of  the  ‘ 

rotations  of  sections  parallel  to  the  direction  of  flight  due  to 
concentrated  unit  torques  in  planes  parallel  to  the  plane  of  symmetry 
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or  concentrated  unit  loads  at  the  reference  line .  These  rotations  in 
radians  are  entered  in  tables  of  the  form: 


M 

WIST  AT  STATION  j/a^  DUE  TO 
UNIT  CONCENTRATED  TORQUE  AT 


\yi /% 
y/sw  \. 

0.2 

0.4 

0.6 

0.8 

0.9 

1.0 

0 

0.2 

0.4 

0.6 

0.8 

0.9 

M 

WIST  AT  STATION  y/sw  DUE  TO 
UNIT  CONCENTRATED  LOAD  AT 


7l/sw 


\yi/% 

yAv 

0.2 

0.4 

0.6 

0.8 

0.9 

1.0 

0 

0.2 

0.4 

0.6 

0.8 

0.9 

These  particular  tables  would  be  used  for  a  six— point  analysis;  similar 
tables  would  be  used  for  a  ten-point  analysis.  In  either  case  it  is 
to  be  noted  that  the  twists  are  measured  at  values  of  y/sw  from  0  to  0.9, 

whereas  the  loads  are  applied  at  yj/s w  values  from  0.2  to  1.0.  The 

tables  obtained  in  this  manner  constitute  the  desired  influence- 
coefficient  matrices. 


If  the  wing  sections  are  found  to  twist  nonuniformly,  so  that  they 
become  cambered  in  effect,  the  angles  of  twist  as  to  be  entered  in 
the  influence-coefficient  matrices  have  to  be  defined  in  a  different 
manner  according  to  whether  the  aeroelastic  analysis  is  performed  for 
subsonic  or  supersonic  speeds.  At  subsonic  speeds  the  lift  depends  on 
the  slope  of  the  mean  camber  line  at  the  three— quarter— chord  point  so 
that  the  effective  angle  of  attack  is  ’ 


1 6 
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At  supersonic  speeds  tlie  lift  depends  primarily  on  the  average  slope  of 
the  mean  camber  line,  so  that 


^LE  ^TE 


<Xa  = 


(31) 


Determination  of  the  Aerodynamic  Parameters 

The  selection  of  the  aerodynamic  parameters  and  e.^  for  the 

calculation  of  the  divergence  speed  has  been  discussed  in  reference  3* 

For  calculating  the  aerodynamic  loading  at  a  given  flight  condition  the 
aerodynamic  parameters  are  chosen  for  that  flight  condition.  The 
effective  lift-curve  slopes  1%  and  nig  are  applicable  only  to  subsonic 

subcritical  speeds.  At  higher  speeds  no  simple  span  correction  is 
available;  neglect  of  the  span  correction  tends  to  be  conservative  for 
calculation  of  the  divergence  speed  and  the  aerodynamic  loading,  however. 


Instructions  for  Solution 


Two  sets  of  integrating  matrices  have  been  prepared,  one  for  a  six- 
point  solution  and  one  for  a  t en — po i nt  solution.  The  former  should  be 
adequate  for  all  practical  purposes;  only  where  the  stiffness  curves  are 
very  irregular  near  the  root  does  the  ten— point  solution  have  to  be 
resorted  to.  The  points  considered  by  the  two  sets  of  tables  are  at 


JL  =  0,  0.2,  0.4,  0.6,  0.8,  and  0.9  for  the  shorter  solution  and 
SA 


-3-=  0,  0.1,  0.2,  0.3,  0.4,  0.5,  0.6,  0.7, 

SA 


0.8,  and  0.9  for  the  longer 


solution.  The  procedure  to  be  followed  for  either  solution  is  identical; 
although  computing  forms  are  presented  in  this  paper  only  for  the  six- 
point  solution,  their  extension  to  apply  to  the  ten-point  solution  is 


obvious. 


Calculation  of  the  matrices.-  The  first  step  in  the  aeroelastic 
analysis  by  means  of  the  stiffness  curves  is  the  calculation  of  the 
aeroelastic  matrix  [A]  from  the  physical  and  geometrical  parameters 
of  the  wing.  These  parameters  are  conveniently  tabulated  in  a  form  of 
the  type  shown  in  table  Yl(a).  The  computation  is  then  carried  out 
according  to  the  instructions  of  table  Vl(b),  each  step  in  the  procedure 
being  identified  by  the  number  in  the  upper  left  comer  of  each  box. 

It  must  be  kept  in  mind  that  many  of  the  operations  call  for  matrix 
multiplications  where  the  order  of  the  multiplicands  is  of  importance. 

(A  "brief  summary  of  matrix  methods  is  presented  in  the  appendix. )  The 
aeroelastic  matrix  is  obtained  as  the  last  step  (step  13 )  of  "the 
computations  in  this  form  which  constitutes  an  evaluation  of  the 
A  matrix  given  in  equation  (15)« 
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A  special  case  arises  when  e-,  is  zero.  If  e,  is  not  zero 


Lr  1 

along  the  remainder  of  the  span,  its  value  at^  some  point  other  than  the 

I 

root  may  be  used  as  a  reference  value.  —  ~ 


The 


1  /_c_ 

’l„  \cr> 


e 

L  tr 


matrix  and  the 


•-  ~  «■>* 

multiplying  factors  of  steps  8  and  9  as  well  as  the  definition  of  the 
parameter  a  are  then  based  on  this  other  reference  value  rather 
than  e^.  If  ej_  is  zero  along  the  entire  span,  step  1  and  steps  3 

to  8  may  be  omitted  and  steps  9  to  13  should  be  modified  as  follows; 


Step  9 


(EE). 


y  jojM  fV  ] 
B*r  A  L  °J 


(GJ)_^  ®w  tan  A 


Step  10 


[®]-  [®] 


Step  11 
Step  12 


As  is 
Omit 


Step  13  Wei=0  -  [©]  [@] 


If  influence  coefficients  of  the  proper  type  are  available,  the 
calculation  of  the  aeroelastic  matrix  lA*]  is  carried  out  directly 
by  means  of  equation  (28). 


Solution  for  divergence  dynamic  pressure.-  In  order  to  determine 
the  value  of  the  parameter  a  or  a1  corresponding  to  divergence,  the 
aeroelastic  matrix  [AJ  or  [A*]  is  iterated  (see  appendix)  as  indicated 
by  equation  (l8).  Table  Tll(d)  may  be  used  for  this  purpose.  The  result 
is  the  critical  value  of  a  or  a*.  The  divergence  dynamic  pressure  is 
then  calculated  from  equation  (l6)  or  (27).  It  is  to  be  noted  that  this 
pressure  will  be  in  pounds  per  square  inch.  Since  the  aeroelastic  matrix 
is  independent  of  the  Mach  number,  except  insofar  as  e^  varies  with  Mach 

number,  the . same  critical  value  of  a  may  be  used  to  calculate  the  diver¬ 
gence  dynamic  pressure  for  an  entire  range  of  Mach  numbers.  If  the  value 
of  e-^  changes,  however,  as  it  does  between  the  subsonic  and  supersonic 

region,  the  critical  value  of  a  has  to  be  calculated  for  both  values 
of.  e^ . 

If  the  value  of  ej  is  zero  along  the  entire  span  and  the  [a] 

matrix  has  been  calculated  according  to  the  modified  instructions, 
iteration  of  the  matrix  will  give  the  value  of  the  parameter  d  at 
divergence.  From  the  definition  of  d  the  divergence  dynamic  pressure 
may  then  be  calculated. 
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Solution  for  aerodynamic  loading.—  In  order  to  calculate  the  aero¬ 
dynamic  loading  corresponding  to  a  given  flight  condition  and  geometric 
angle— of-attack  distribution  the  aeroelastic  matrix  [A]  or  [A*]  .  is 
multiplied  hy  the  value  of  a  or  a*  calculated  for  the  given  flight 
condition  and  subtracted  from  the  unit  matrix  [IJ*  (See  equation  (17)°) 
The  result  may  be  entered  in  table  Yll(b).  Again  it  must  be  noted  that  the 
value  of  the  aeroelastic  matrix  varies  with  the  flight  condition  if  e-^ 
varies,  so  that  the  aeroelastic  matrix  corresponding  to  the  proper  ex 
value  must  be  selected.  The  resulting  matrix  constitutes  the  coefficients 
of  a  set  of  simultaneous  linear  algebraic  equations  for  the  unknown  values 

of  the  effective  angle— of— attack  distribution  of  the  deformed  wing 

>  in  terms  of  the  known  angle— of— attack  values  of  the  rigid 
L6!  J 

wing  |ag|.  Table  Yll(b)  is  set  up  for  the  calculation  of  the  additional 

loading,  the  damping-in-roll  loading,  and  a  third  arbitrary  loading;  as 
many  loadings  as  desired  may,  of  course,  be  calculated  by  this  method. 

The  solution  of  the  equations  may  be  carried  out  in  any  convenient  manner. 
The  form  of  table  Vll(b)  has  been  prepared  for  use  in  conjunction  with 
Grout* s  method  of  solving  linear  simultaneous  equations  (reference  6). 

In  the  case  where  ej  is  zero  along  the  span,  the  headings  at  the 
top  of  table  Vll(b)  Bhould  be  modified  to  read 


[p]  -  *Mei=0] 


■wiier©  [A]e  has  "been  calculated  according  to  the  modified  instructions 
and  d  has  been  obtained  by  iterating  [A]e^_Q. 


The  values  of 


&»} 


calculated  for  the  additional  load 


distribution*  (ag  =  l)  constitute  values  of  the  ratio  °l^/clrw 
or  (cci)^  cci)  in  view  of  the  assumptions  made  concerning  the  air 

forces.  The  section  loading  of  the  flexible  wing  is  obtained  from  the 


relation 


cci =  c^iU!rE 


(32) 
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a 


or  in  dimensionless  form 


221  =  m 

Cr  G± 


(32a) 


The  wing  lift  coefficient  defined  by  the  relation 

<33) 

and  the  wing  tending-moment  coefficient  defined  by 

may  be  obtained  by  integrating  the  load  distribution.  These  integrations 
may  be  performed  conveniently  by  multiplying  the  cc2/cr  values  by  the 

first  rows  of  the  |KqJ  and  matrices,  respectively.  Thus 


and 


(35) 


(35a) 


r  -  8wcr  Bw 


(56) 
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y 

The  lateral  center  of  pressure  of  the  wing  load  —  may  then  he 

sw  . 

determined  from  the  relation 


7_  _  21b  CbMw 
BW  8w  °Iv 


(37) 


The  fore-and-aft  location  of  the  aerodynamic  center  of  the  wing  load 
measured  rearward  of  the  leading  edge  of  the  mean  aerodynamic  chord  as 
a  fraction  of  the  mean  aerodynamic  chord  may  he  estimated  from  the 
relation 


(ac)w 

CMAC 


ac  + 


7  7MAC 
cMAC 


tan  Aac 


(38) 


where  A  is  the  sweep  of  the  section-aerodynaml c~center  line, 
ac 

For  any  other  geometrical  angle-of-attack  distributions  such  as 
those  due  to  built— in  twist  or  those  due  to  rolling,  the  same  section 
lift-curve  slope  should  he  used  as  for  the  structural  deformations,  so 

that  mg  is  replaced  by  nig  and  — •  is  unity  in  equations  (32),  (35) 

and  (36).  For  the  damping-in-roll  distribution  with  a  tip  helix  angle 
of  1  radian 


(39) 


The  rolling-moment  coefficient  due  to  the  wing  load  is  defined  by 


C 


-  El 

hr  qSb 


(40) 


It  is  seen  to  be  twice  the  wing  bending-moment  coefficient. 

The  contribution  of  the  wing  to  other  stability  derivatives  may 
be  obtained  similarly  by  integrating  the  load  distributions  due  to  the 
angle-of— attack  distributions  caused  by  the  motion  under  consideration, 
as  described  in  reference  7;  in  the  case  of  swept  wings,  particular  care 
must  be  taken  in  selecting  the  proper  angle-of-attack  distribution  and 
in  accounting  for  the  lateral  inclination  of  the  lift  vector.  (See 
reference  4.) 

If  the  aerodynamic  loading  or  the  stability  derivatives  are  to  be 
obtained  for  a  wide  variety  of  flight  conditions,  it  is  convenient  to 
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systematize  the  calculations  in  the  following  manner;  The  aeroelastic 

matrix  is  computed  for  hoth  the  subsonic  and  supersonic  aerodynamic _ 

center  values  and  iterated  for  both  cases  to  obtain  the  subsonic  and 
supersonic  values  of  the  divergence  parameter  ap.  Prom  these  values 

the  divergence  dynamic  pressure  may  be  computed  by  means  of  equation  (15) 
and  plotted  against  Mach  number,  as  suggested  in  reference  3;  on  the 
same  plot  values  of  the  actual  dynamic  pressure  may  be  plotted  against 
Mach  number  for  various  altitudes  of  interest.  Such  a  plot  for  a  wing, 
the  physical  characteristics  of  which  are  given  in  figure  3,  is  shown  in 
figure  4. 

Since  at  a  given  Mach  number  the  ratio  a/ap  is  equal  to  the 
ratio  q/qp,  the  range  of  a/a^  values  of  interest  may  be  established 

from  this  plot  for  both  the  subsonic  and  the  supersonic  region.  Several 
representative  a/a^  values  may  then  be  chosen  within  the  given  ranges 
and  the  corresponding  values  of  a  computed  from  the  previously  calculated 
ap  values.  The  aerodynamic  loading  is  calculated  for  these  values  of  a 
using  the  appropriate  [a]  matrix  and  plotted  in  the  form  of 
(cc2 ^fv/(CC2  rw*  ra^^-°  a/ap  as  a  parameter.  Prom  these 

curves  or  from  the  jg  a  values  the  lift  coefficients  may  be  obtained 

®1 

cL^rw  821(1  plotted  against  a/aD  or  q/qD;  the 
other  coefficients  may  be  obtained  and  plotted  in  a  similar  form. 

For  any  specific  flight  condition  the  value  of  a/aD  may  then 
be  obtained  from  the  plot  of  q  and  qp  against  Mach  number.  The 

loading,  lift  coefficient,  or  other  item  of  interest  may  be  obtained 
from  the  plots  which  give  these  items  in  terms  of  the  rigi dewing  values. 
Once  the  rigid-wing  values  at  the  given  Mach  number  are  known,  the 
flexible-wing  values  may  then  be  obtained  immediately. 


Illustrative  Example 

In  order  to  illustrate  the  method  described  in  the  preceding 
sections,  a  typical  swept  wing  has  been  analyzed.  The  physical  and 
geometrical  parameters  of  the  wing  are  shown  in  figure  3  and  the  upper 
part  of  table  YIH  (which  follows  the  form  of  table  VT(a)).  The 
chord,  eqc  ,  and  stiffness  matrices  have  been  obtained  from  the 
given  parameters  and  are  shown  in  the  lower  part  of  table  VIII. 

The  calculation  of  the  aeroelastic  matrix  for  the  subsonic  case 
has  been  carried  out  by  means  of  the  form  of  table  VT(b).  All  but 
three  of  the  steps  of  the  computation  are  shown  in  table  IX  numbered 
in  the  same  order  as  in  table  Vl(b).  Steps  1,  2,  6,  7,  11,  and  12 
constitute  matrix  multiplications,  which  are  carried  out  in  the  order 
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indicated;  steps  5  and  13  constitute  matrix  additions  or  subtractions; 
steps  3  and  k  constitute  multiplications  of  matrices  "by  constants. 

Tha  aeroelastic  matrix  is  iterated  in  table  X(a)  (which,  follows 
the  form  of  table  Yll(a))  to  yield  a  value  of  a^  =  -2.208.  Prom  this 

value  and  a  value  of  computed  in  the  same  manner  for  supersonic 

speeds,  the  divergence  dynamic  pressure  has  been  calculated  by  means  of 
equation  (l6)  on  the  basis  of  estimated  values  of  the  effective  lift- 
curve  slope.  The  variation  with  Mach  number  of  the  divergence  dynamic 
pressure,  the  actual  dynamic  pressure  at  sea  level,  and  the  estimated 
effective  lift-curve  slope  is  shown  in  figure  4. 

For  a  value  of  —  =  -0.25,  such  as  would  be  obtained  approximately 
aD 

at  a  Mach  number  of  1.0,  the  aerodynamic  loading  has  been  calculated  for 
the  additional-angle-of-attack  case  and  the  damping-in— roll  case  in 
table  X(b),  which  follows  the  form  of  table  VH(b).  The  values  of  ag 

for  the  damping— in— roll  case  have  been  calculated  from  equation  ( 39 ) • 

The  aerodynamic  loadings,  in  addition  to  those  calculated  for  other 
q/q^  values,  have  been  plotted  in  figure  5  as  ratios  of  the  flexible- 
wing  loadings  to  the  rigid— wing  loadings.  The  curves  have  been  integrated 
to  yield  wing  lift  and  rolling-moment  coefficients  as  well  as 
the  aerodynamic  center  of  the  wing  load,  which  are  shown  in 

table  X'(b )  for  the  case  of  — —  =  —0 . 25  and  which  are  plotted  against 

aD 

— —  in  figure  6. 

ID 

The  wing  lift  coefficient  is  defined  in  such  a  manner  that  if 
the  fuselage  lift  is  known  and  made  dimensionless  by  dividing  by  q 
and  S  the  resulting  fuselage  lift  coefficient  may  be  added  directly 
to  the  wing  lift  coefficient.  This  definition  and  the  fact  that 
figure  5(a)  is  plotted  over  the  fraction  of  the  wing-alone  span  sw 

explain  the  fact  that  the  area  under  the  curve  of  figure  5(a)  is  not  1. 
Th9  aerodynamic  center  as  plotted  in  figure  6  constitutes  the  center 
of  pressure  of  only  the  wing  load.  In  order  to  obtain  the  airplane 
aerodynamic  center,  the  magnitude  and  center  of  pressure  of  the  fuselage 
load  would  have  to  be  known  and  taken  into  account. 


DISCUSSION 


Both  the  aerodynamic  and  the  structural  assumptions  made  in  this 
analysis  are  more  realistic  than  those  made  in  reference  3 ■  The  device 
employed  in  this  analysis  of  calculating  the  air  forces  for  wing  sections 
parallel  to  the  direction  of  flight  and  then  transferring  them  to 
sections  perpendicular  to  the  elastic  axis  obviates  the  necessity  of 
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replacing  the  actual  wing  with  one  the  root  and  tip  of  which  are 
perpendicular  to  the  elastic  axis  for  the  purpose  of  analysis.  Further¬ 
more,  the  inclusion  of  Falloierts  modification  (see  appendix)  in  the 
integrating  matrices  has  the  effect  of  rounding  off  the  load  distribu¬ 
tion  approximately  in  the  manner  observed  at  subsonic  speeds.  At 
supersonic  speeds  the  load  distributions  do  not  go  to  zero  in  the  manner 
assumed  in  Fa1Tmerts  modification,  but  even  at  supersonic  speeds  there 
is  some  reduction  of  load  at  the  tip,  the  total  magnitude  of  which  is 
not  far  from  the  reduction  obtained  by  Falkner’s  modification. 


Only  one  aerodynamic  assumption  is  still  made:  that  induction 
effects  may  be  approximated  by  an  over— all  reduction  of  the  strip  theory 
loading  (rounded  off  as  previously  described)  at  subcritical  speeds  and 
may  be  neglected  at  supersonic  speeds.  The  effects  of  aerodynamic 
induction  could  be  taken  into  account  more  accurately  by  using  aero¬ 


dynamic  influence— coefficient  matrices  instead  of  the  effective  lift- 


curve— slope  concept  and  the 


0 

°  21 

K 

j  and 

matrices  used  in 


this  analysis.  Available  methods  of  calculating  such  influence  coeffi¬ 
cients  for  wings  of  arbitrary  plan  form  at  subsonic  and  supersonic  speeds 
particularly  those  suitable  for  wings  with  large  amounts  of  sweep,  are 
either  too  inaccurate  or  too  time-consuming  for  practical  purposes, 
however . 


Although  the  analysis  of  this  paper  has  been  performed  for  wings 
consisting  of  uncambered  sections,  it  is  directly  applicable  as  well  to 
the  determination  of  the  additional  loading  of  wings  with  cambered 
sections.  The  loading  of  such  wings  due  to  the  section  pitching  moment 
at  zero  lift  may  be  determined  by  modifying  the  analysis  somewhat. 

The  assumption  of  an  effective  root  perpendicular  to  the  elastic 
axis  made  in  reference  3  for  the  purposes  of  calculating  the  structural 
response  is  carried  over  in  this  analysis.  It  is  modified,  however, 
to  the  extent  that  the  root  is  no  longer  considered  to  be  rigid  as  in 
reference  3,  but  flexible,  both  in  torsion  and  bending.  It  has  been 
demonstrated  in  reference  5  that  the  deflections  of  a  swept  beam  may  be 
estimated  on  that  assumption,  provided  the  root— rotation  parameters  are 
known.  By  assuming  the  effective  root  at  the  intersection  of  the  elastic 
axis  with  the  side  of  the  fuselage,  the  root  bending  due  to  bending 
moment  and  root  twist  due  to  torque  are  minimized.  The  bending  due  to 
twist  and  twist  due  to  bending  are  the  same  regardless  of  the  location  of 
the  effective  root. 

The  method  of  Introducing  the  root  rotations  into  the  analysis  by 
means  of  the  [iqj  matrix  assures  that  the  structural  twist  in  planes 
parallel  to  the  direction  of  flight  is  zero  at  the  fuselage.  From 
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figure  2  it  is  seen  tlifct  the  local  values  of  the  root  rotation  either 
tend  to  approach  zero  at  the  root  or  tend  to  cancel  each  other.  If  the 
root — rotation  constants  are  known,  the  structural  deforxnations  can  there 
fore  he  predicted  quite  accurately  hy  the  assumptions  made. 

The  manner  in  which  the  equations  of  equilibrium  are  solved  hy 
means  of  the  integrating  matrices  accounts  for  the  true  chord  and 
stiffness  variations.  It  does  not  necessitate  replacement  of  the 
actual  wing  hy  constant-chord  segments  with  all  the  flexibility 
concentrated  at  the  ends  of  the  segments,  an  approach  which  has  been 
used  extensively  in  the  work  on  aeroelastic  problems  of  straight  wings. 

A  further  refinement  which  obviates  the  necessity  for  making  any 
structural  assumptions  other  than  that  of  small  deflections  is  the  use 
of  measured  influence  coefficients  in  the  aeroelastic  analysis.  Wherever 
such  coefficients  are  available  it  is,  of  course,  of  advantage  to  use  them. 

No  explicit  account  has  been  taken  in  the  analysis  of  the  effects 
of  the  inertia  loading  on  the  structural  deformations  and  hence  the 
aerodynamic  loading.  On  swept  wings,  in  particular,  their  effects 
may  be  considerable.  For  the  purposes  of  this  analysis,  however,  the 
structural  deformations  due  to  inertia  loading  may  be  considered  part 
of  the  geometric  angle  of  attack  and  the  rigid-wing  geometric  angle  of 
attack  may  be  modified  accordingly.  The  deformations  due  to  the  inertia 
loading  may,  incidentally,  be  calculated  conveniently  by  means  of  the  Kq, 
Kg,  and  Kg  matrices. 

Some  of  the  general  observations  made  in  reference  3  concerning  the 
divergence  phenomenon  are  corroborated  by  the  example.  As  expected  of 
a  wing  with  a  considerable  amount  of  sweepback,  the  divergence  dynamic 
pressure  is  negative.  Consequently  the  wing  cannot  diverge.  The 
divergence  dynamic  pressure  is  useful  as  a  reference  value,  however; 
the  values  of  the  load  distribution  and  the  stability  parameters  divided 
either  by  the  corresponding  rigid-wing  values  or  by  the  section  lift- 
curve  slope  depend  only  on  the  ratio  of  the  actual  to  the  divergence 
dynamic  pressure. 

The  type  of  plot  shown  in  figure  4  is  therefore  quite  useful  in 
the  analysis  of  aeroelastic  phenomena.  As  pointed  out  in  reference  3, 
this  chart  may  also  be  used  to  estimate  the  actual  divergence  dynamic 
pressure  where  there  is  a  possibility  that  the  wing  may  diverge.  It 
appears  that  the  critical  values  will  occur  at  either  extremity  of  the 
transonic  region.  In  the  transonic  region  proper  the  lift-curve  slope 
usually  appears  to  be  lower  than  at  the  extremities ,  so  much  so  that  the 
decrease  in  lift— curve  slope  even  tends  to  overbalance  any  forward 
shift  in  aerodynamic  center. 

As  would  be  expected  qualitatively,  the  effect  of  wing  flexibility 
in  the  case  of  the  example  wing  is  to  unload  the  wing  tips  owing  to  the 
fact  that  they  bend  up.  The  lift  carried  by  the  wing  is  therefore  less 
than  that  carried  by  a  rigid  wing,  the  center  of  pressure  being  farther 
inboard  and  the  aerodynamic  center  farther  forward. 
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®ie  difference  between  the  supersonic  and  subsonic  values  of  the 
.  loading,  the  lift  and  rolling-moment  coefficients,  and  the  aerodynamic 
center  for  a  given  value  of  a/a^  is  due  to  the  difference  in  the 

e2  distributions.  If  the  distributions  were  the  same,  the  subsonic 
and  supersonic  variations  of  these  quantities  with  a/aD  would  coincide 
despite  the  difference  in  the  e-j_  values. 


Another  item  of  possible  interest  is  the  fact  that  the  variations 
of  the  parameters  dp  and  ap  for  the  example  problem  are  approximately 

linear  (see  fig.  7),  as  would  be  expected  from  the  results  of  the 
analysis  of  reference  3 •  The  deviations  from  linearity  are  most  pronounced 
near  the  points  for  d  =  0  (that  is,  A  =  0°).  They  are  due  to  the  effects 
of  the  root  rotations,  in  particular,  the  bending  due  to  torsion  and  torsion 
due  to  bending;  these  effects  were  neglected  in  the  approximate  analysis  of 
reference  3 •  The  points  of  figure  7  correspond  to  the  example  wing  and 
the  wings  that  would  be  obtained  by  rotating  the  example  wing  to  the 
unswept  and  37*5°  sweptforward  positions  in  such  a  manner  as  to  keep  the 

r*  />  ao  A  / _ \  * 


elrc£  008  A  (El) 
s, 


as  well  as  the  chord,  stiffness,  and 


parameters  — .  - 

8a  «sr>P 

moment-arm  (ex)  distributions  constant.  Points  are  shown  for  both  the 
subsonic  and  supersonic  variations  as  well  as  for  the  case  when  e-,  =  0 
over  the  entire  span  (aD  =  0).  The  difference  between  the  subsonic  and 


supersonic  lines  is  due  entirely  to  the  difference  in-  the  ej_  distribu¬ 
tion;  if  the  distributions  were  the  same,  as  would  be  the  case  if  the 
elastic  axis  were  at  a  constant  fraction  of  the  chord,  the  variations 
would  be  the  same,  regardless  of  the  difference  in  the  e1  values. 


The  present  analysis  is  concerned  only  with  wing  or  tail  loads; 
the  total  loads  are  obtained  by  adding  the  fuselage  loads  (which  may 
be  assumed  to  be  unaffected  by  flexibility)  to  the  wing  or  tail  loads 
obtained  from  the  analysis.  The  amount  of  load  carried  by  a  flexible 
wing  and  the  manner  of  its  distribution  can  consequently  be  estimated  by 
the  method  presented  herein  if  the  contribution  of  the  fuselage  is  known 
at  low  dynamic  pressures,  that  is,  for  the  "rigid-^wing"  case. 


The  fuselage  has  a  considerable  effect  on  some  of  the  stability 
parameters  as  well,  although  in  the  case  of  others,  such  as  C ^  ,  the 

effect  is  negligible.  Other  effects  that  may  have  to  be  accounted  for 
in  calculating  stability  derivatives  are  the  boundary— layer  behavior 
and  tip  suction.  The  boundary— layer  effect  may  be  accounted  for  by 
using  a  section  lift-curve  slope  corrected  for  boundary— layer  effects 
to  calculate  the  angle— of-attack  distribution  of  the  flexible  wing  at 
the  flight  conditions  of  interest  and  then  obtaining  the  lift  and  drag 
distributions  corresponding  to  that  angle-of-attack  distribution. 
Lateral  tip  suction  may  be  important  on  low— aspect -ratio  and  highly 
swept  wings.  Since  it  does  not  affect  the  lift  distribution,  it  may  be 
taken  into  account  by  calculating  the  angle-of-attack  distribution  of 
the  flexible  wing  and  estimating  the  tip  suction  corresponding  to  the 
actual  angle  of  attack  at  the  tip. 
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In  calculating  stability  derivatives  it  is  well  to  keep  in  mind 
that  the  method  presented  in  this  paper  is  based  on  a  modified  strip 
theory,  unless  aerodynamic  influence-coefficient  matrices  are  used.  The 
calculated  derivatives  may  therefore  be  somewhat  in  error,  particularly 
if  in  calculating  them  the  moment  of  a  load  distribution  has  to  be 
determined.  If  there  is  reason  to  suspect  that  the  modified  strip 
theory  is  inadequate  for  calculating  a  giten  derivative,  the  derivative 
may  be  calculated  for  the  rigid-wing  case  by  a  more  refined  method;  the 
results  calculated  by  the  method  of  this  paper  may  then  be  used  to 
correct  the  accurate  rigid-wing  value  for  the  effect  of  structural 
flexibility. 


CONCLUDING  REMAKES 


A  method  has  been  presented  for  calculating  the  aerodynamic  loading, 
the  divergence  speed,  and  certain  stability  derivatives  of  wings  and 
tail  surfaces  of  arbitrary  plan  form  and  stiffness.  Provisions  have 
been  made  for  using  either  stiffness  curves  and  root-rotation  constants 
or  influence  coefficients  in  the  structural  part  of  the  analysis.  Strip 
theory  with  over-all  reduction  and  rounding  off  at  the  tip  to  take 
account  of  aerodynamic  induction  have  been  used  for  the  aerodynamic  part 
of  the  analysis.  Computing  forms,  tables  of  numerical  constants  required 
in  the  analysis,  and  an  illustrative  example  are  included  to  facilitate 
calculations  by  means  of  the  method. 
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APPENDIX 

* 

SUMMARY  OP  MATRIX  ALGEBRA  PERTINENT  TO  THE  ANALYSIS 


For  the  convenience  of  the  reader  unfamiliar  with  matrix  terminology 
a  summary  of  matrix  definitions  and  methods  is  presented  in  the  following 
sections.  For  a  more  complete  discussion  of  matrix  methods  the  reader 
is  referred  to  any  text  on  matrices,  for  instance  reference  8. 


Definitions 

A  matrix  is  a  rectangular  array  of  numbers,  called  elements,  written 
down  in  rows  and  columns.  A  column  matrix  consists  of  a  single  column , 
a  row  matrix  of  a  single  row.  A  square  matrix  has  as  many  rows  as  it 
has  columns.  The  diagonal  of  a  square  matrix  from  the  upper  left  to  the 
lower  right  is  called  the  principal  diagonal.  A  matrix  all  the  elements 
of  which  are  zero  except  for  those  on  the  principal  diagonal  is  called  a 
diagonal  matrix.  If  all  of  these  elements  are  unity,  the  matrix  is 
termed  the  unit  matrix. 


Matrix  Algebra 

Addition.—  Two  matrices  can  be  added  or  subtracted  if  both  have  the 
same  number  of  rows  and  columns.  The  addition  or  subtraction  is  carried 
out  by  adding  to  or  subtracting  from  each  element  of  the  first  matrix 
the  corresponding  element  of  the  second  matrix. 

Multiplication  by  a  constant.-  A  matrix  is  multiplied  by  a  constant 
by  multiplying  each  element  by  that  constant. 

Matrix  multiplication.-  Two  matrices  can  be  multiplied  by  each  other 
if  the  second  has  as  many  rows  as  the  first  has  columns.  The  elements 
of  the  resulting  matrix  are  obtained  by  multiplying  the  elements  in  the 
corresponding  row  of  the  first  matrix  by  those  of  the  corresponding 
column  of  the  second  matrix  in  the  following  order:  The  first  element 
of  the  row  is  multiplied  by  the  first  element  of  the  column,  the  second 
by  the  second,  and  so  forth.  The  Bum  of  the  products  obtained  in  this 
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manner1  is  the  value  of  the  element  of  the  product  matrix.  Schematically 
this  process  may  he  illustrated  as  follows : 


Q=aA  +  bB  +  cC  +  cLD  +  eE  +  fF 


It  must  he  emphasized  that  in  multiplying  matrices  by  each  other  their 
order  is  of  importance.  As  the  two  matrices  under  consideration  are 
written  the  matrix  at  the  left  (the  m  matrix)  is  said  to  he  post- 
multiplied  hy  the  other,  (the  M  matrix);  or  the  M  matrix  may  he  said 
to  he  premultiplied  hy  the  m  matrix,  in  order  to  distinguish  the  manner 
in  which  they  are  multiplied.  If  the  two  matrices  were  written  in  the 
reverse  order  and  then  multiplied  according  to  the  foregoing  instructions, 
that  is,  if  the  [m]  matrix  were  postmult iplied  hy  the  [mj  matrix,  the 
element  of  the  second  row  and  fourth  column  of  the  product  matrix  |MJLmJ 
would  clearly  not  have  the  vfilue  Q  in  general,  nor  would,  in  general, 
any  other  element  have  the  value  it  would  have  if  the  two  matrices  were 
multiplied  in  the  order  shown.  Consequently  it  is  important  to  observe 
the  order  in  which  the  matrices  are  written  down  in  the  computing 
instructions. 

Matrix  iteration.—  The  purpose  of  iterating  a  square  matrix  is  to 
determine  the  column  matrix  or  matrices  which,  if  postmultiplied  hy  the 
given  square  matrix,  yield  the  same  column  matrix  except  for  a  constant 
multiplier.  It  is  the  value  or  values  of  these  multipliers  which 
constitute  the  desired  characteristic  values  of  the  matrix. 

The  iteration  is  carried  out  hy  assuming  a  "trial  column  (the 
column  shown  in  table  VII  is  convenient  for  the  purpose  of  this  analysis) 
•and  premultiplying  it  hy  the  given  square  matrix  to  yield  a  "result" 
column.  Thft  elements  of  the  result  column  including  the  last  are 
divided  hy  the  last  element  of  the  result  column  and  entered  as  a 
second  trial  column.  The  second  trial  column  is  then  premultiplied 
hy  the  square  matrix  to  yield  a  second  result  column.  The  procedure  is 
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repeated  until  the  same  value  (within  the  desired  accuracy)  is  obtained 
twice  in  succession  for  the  last  element  of  the  result  matrix.  The 
reciprocal  of  this  value  is  the  desired  (lowest)  characteristic  value 
of  the  matrix,  that  is,  the  lowest  critical  value  of  a^,  in  the 
analysis  of  this  paper. 

Another  way  of  estimating  a  first  trial  column  in  this  analysis 
is  to  add  the  elements  in  each  row  of  the  A  matrix,  enter  the  six  sums 
in  the  first  result  column,  and  treat  them  as  if  they  had  "been  obtained  by 
multiplying  the  A  matrix  by  a  first  trial  column. 


Derivation  of  the  Integrating  Matrices 

Although  familiarity  with  the  derivation  of  the  integrating 
matrices  is  not  essential  to  the  application  of  the  method  of  this 
paper,  an  outline  of  the  derivation  is  presented  because  of  its  general 
interest . 

The  integrating  matrices  used  in  this  paper  are  based  on  the  same 
concept  as  Simpson* s  rule  —  replacement  of  the  actual  function  which  is 
to  be  integrated  by  parabolic  segments.  If  the  function  y  has  the 
values  yn_2_,  yn,  and  yn+p,  respectively,  at  the  equally  spaced 

points  x^,  xn,  and  xn+i,  the  following  relations  are  seen  to  be 

true  for  a  second— degree  parabola  passed  through  the  three  known  points: 


y  =  Fn  +  |(yn+l  -  y-n-lX*  -  xn)  +  |<7n+l  “  2?n  +  Fn-lX*  “  (Al) 
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J  (x  -  s^)y  dx  =  (J-  ^  Ax^y^  +  (°)Jn  +  (3  ^^n+l  (A5^ 

xn-l 


fn+1  (x  -  *n>  dx  =  ^  7n_i  +  (|  Ax^)yn-1  +  Ql  ^X+l  (A6) 

% 


where 


Ax  =  -  Xq_i  =  %+x  ~  % 

The  different  integrations  over  the  parabolic  segments  may  thus 
he  performed  by  multiplying  the  given  y  values  by  the  multiplying 
factors  indicated  in  equations  (A2)  to  (A 6). 

Since  load  distributions  at  subsonic  speeds  go  to  zero  with  infinite 
slope  at  the  tip  and  the  ordinary  second-degree  parabola  furnishes  a 
poor  approximation  to  such  a  distribution,  V.  M.  Falkner  has  suggested 
that  a  curve  of  the  type 

y  =  Aq  +  A1(l  -  x)^2  +  Agd  -  i)3/2  (A0) 

be  passed  through  the  last  three  points  of  the  load-distribution  curve 
at  the  tip  (x  =  l).  On  the  basis  of  the  approximation,  relations 
equivalent  to  equations  (Al)  to  (A 6)  may  be  derived.  The  multiplying 
factors  for  the  last  two  segments  are  then  based  on  these  equivalent 
expressions  rather  than  those  of  equations  (A2)  to  (a6). 

The  integrating  factors  of  equations  (A2)  to  (A 6)  may  be  assembled 
directly  into  integrating  matrices.  The  Kg  matrix,  for  instance,  is 

set  up  to  perform,  the  integration  JXj  dx.  If  at  the  upper  limit  x  =  0.1 

and  if  the  ten-point  matrix  (table  111(b))  is  to  be  used,  the  factors  0.04167, 
0.06667,  and  -O.OO833  may  be  obtained  from  equation  (A4)  since  x^  -  0, 
xn  =  0.1,  and  Ax  =  0.1;  similarly,  if  for  the  same  case  the  integration 
is  extended  to  xn+i  =  0.2  as  the  upper  limit,  the  integrating  factors 
0.03333,  0.13333,,  and  0.03333  will  be  obtained  from  equation  (A2). 

These  factors  coiistitutb ;  the  second  and  third  rows  of  the  matrix  K3; 
since  the  integrations  are  independent  of  the  y  values  other  than  the 
first  three/  the  other  y  'values  are  multiplied  by.  zero  in  these  two 
rows.  In  order  to  extend  the  integration  to  x  =  0.3  an  integration  is 
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again  performed  up  to  x  =  0.2  and  another  integration,  using  another 
parabolic  segment,  is  performed  from  x  =  0.2  to  x  =  0.3.  For  the 
latter  integration  x^  =  0.2,  3^  =  0.3,  and  Ax  =  0.1,  so  that 
equation  (Ah)  again  yields  the  factors  0.04167,  O.O6667,  and  -O.OO833. 
The  y  value  at  x  =  0.2  is  therefore  assigned  a  multiplying  factor 
of  O.O3333  by  the  first  integration  and  a  factor  of  0.04l67  by  the 
second,  or  a  total  factor  of  0.07500.  The  resulting  factors  are 
entered  in  the  fourth  row  of  the  matrix.  All  other  rows  are 
obtained  in  a  similar  manner. 


The  Kj  matrix  is  set  up  to  perform  the  integration  /  y  dx.  The 

Jx 

values  of  the  last  row  of  the  ten— point  Z^  matrix  (table  1(h))  are 

obtained  from  Falimer’s  equivalent  of  equation  (A3)  for  the  curves 
assumed  in  reference  5,  with  x^  =  0.8,  x&  =  0.9,  :%+1  =  1.0,  and 

Ax  =  0.1.  Only  the  multiplying  factors  for  the  y  values  at  x  =  0.8 
and  x  =  0.9  are  listed,  since  the  y  value  at  x  =  1.0  (the  wing  tip) 
is  assumed  to  be  zero  in  this  analysis,  so  that  its  multiplying  factor 
is  immaterial.  The  values  of  the  last  row  but  one  are  obtained  similarly 
from  Falkner,s  equivalent  of  equation  (A2).  The  values  of  the  row 

for  =  0.7  are  obtained  by  using  equation  (A3)  in  the  interval 
“A 

x  =  0.6  to  x=0.8  and  Falkner*s  equivalent  of  equation  (A2)  in  the 

interval  x  =  0.8  to  1.0.  Similarly  the  row  for  -3-  =  0.6  is 

SA 

obtained  by  combining  the  results  of  equation  (A2)  for  the  interval  x  =  0.6 
to  0.8  with  FaHner* s  equivalent  of  equation  (A2)  for  the 
interval  x  =  0.8  to  1.0.  All  other  rows  are  obtained  in  a  similar 
manner . 


Tbe  Kg  matrix  is  set  up  to  perform  the  integration  J  (x  —  xQ)y  dx, 

x0 

where  x  is  the  variable  of  integration  and  x  the  value  of  x  at 
the  lower  limit.  In  applying  the  integrating  factors  of  equations  (A 2) 
through  (A6)  to  this  integration  it  must  be  realized  that 


x0)y  dx  =  (xn  -  x0)  /y  dx  +  /(x  -  xjy  dx 


so  that  the  integrating  factors  for  this  integration  would  be  obtained 

by  adding  (xn  -  xQ)  times  the  factors  of  equation  (A2)  or  (A3)  to 

the  factors  of  equation  (A5)  or  (A6),  respectively,  depending  on  the 

limits  of  the  integration.  The  factors  for  the  different 

segments  (x  =  0.8  to  1.0,  0.6  to  0.8,  and  so  forth)  are  then  combined 

for  any  given  row  (with  its  given  value  of  x0)  in  the  manner  indicated  for 

for  the  Ki  matrix  to  yield  the  K2  matrix. 
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The  matrix  sums  up  the  torques  outboard  of  a  given  point, 
while  the  matrix  gives  the  sum  of  the  moments  of  forces  applied 

outboard  of  a  given  point.  Neither  requires  any  integrations  in  the 
sense  of  equations  (A2)  to  (a6).  For  the  six— point  method  these  two 
matrices  are: 


E4 


*5 


y/sw 

0.2 

0.4 

0.6 

0.8 

0.9  1.0 

0 

0.2 

0.4 

0.6 

0.8 

0.9  1.0 

0.2 

0 

.2 

.4 

.6 

.7  .8 

0.4 

0 

0 

.2 

.4 

.5  -6 

0.6 

0 

0 

0 

.2 

.3  .4 

CO 

• 

O 

0 

0 

0 

0  , 

.1  .2 

0.9 

0 

0 

0 

0 

0  .1 

It  will  be  noted  that  the  moment  arms  which  comprise  the  El  matrix  are 
fractions  of  sv,  so  that  the  matrix  must  be  multiplied  by  the  length  av 

in  order  to  yield  actual  moments  as  stated  in  equation  (21 ). 
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TABLE  I.-  VALUES  OF  THE  INTEGRATING  MATRIX 


(a)  Six— Point  Solution 


n/sA 

0 

.a 

.4 

.6 

.8 

.9 

0 

0.06667 

0.26667 

0.13333 

0.26667 

0.09333 

0.15085 

m2 

-.01667 

.13333 

.15000 

.26667 

.09333 

.15085 

A 

0 

0 

.06667 

.26667 

.09333 

.15085 

.6 

0 

0 

-.01667 

.13333 

.11000 

.15085 

.8 

0 

0 

0 

0 

.02667 

.15085 

•9 

0 

0 

0 

0 

-.01886 

.09333 

(t )  Ten— Point  Solution 


l/sA 

0 

.1 

.2 

.3 

■1 

.5 

.6 

H 

.8 

.9 

0 

0.03333 

0.06667 

0.13333 

0.06667 

0.13333 

0.06667 

0.13333 

0.06000 

0.15085 

.1 

-.00833 

.06667 

.07500 

.13333 

.06667 

.13333 

.06667 

.13333 

.06000 

.15085 

.2 

0 

0 

.03333 

.13333 

.06667 

.13333 

.06667 

.13333 

.06000 

.15085 

.3 

0 

0 

-.00833 

.06667 

.07500 

.13333 

.06667 

.13333 

.06000 

.15085 

m 

0 

0 

0 

0 

.03333 

.13333 

.06667 

.13333 

.06000 

.15085 

.5 

0 

0 

0 

0 

-.00833 

.06667 

.07500 

.13333, 

.06000 

.15085 

.6 

0 

0 

0 

0 

0 

0 

.03333 

.13333 

.06000 

.15085 

.7 

0 

0 

0 

0 

0 

0 

-.00833 

.06667 

.06833 

.15085 

.8 

0 

0 

0 

0 

0 

0 

0 

0 

.02667 

.15085 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

-.01886 

•09333 
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TABLE  II.-  VALUES  OF  THE  INTEGRATING-  MATRIX 


py 


(a)  Six— Point  Solution 


^A 

0 

.2 

.4 

.6 

'  .8 

.9 

0 

0 

0.05333 

0.05333 

0.16000 

0.07314 

0.13792 

.2 

-.00167 

.01000 

.02500 

.10667 

.05448 

.10775 

.4 

0 

0 

0 

•05333 

.03581 

.07758 

.6 

0 

0 

-.00167 

.01000 

.01548 

.04741 

.8 

0 

0 

0 

0 

-.00152 

.01724 

•9 

0 

0 

0 

0 

-.00108 

.00419 

(b)  Ten-Point  Solution 


i/ba 

0 

.1 

.2 

.3 

.4 

.5 

,6 

.7 

.8 

.9 

a 

0.013333 

0.013333 

0.040000 

0.026667 

0.066667 

0.040000 

0.093333 

0.046476 

0.137920 

D 

-.000417 

.002500 

.006251 

.026667 

.020000 

.053333 

.033333 

.080000 

.040476 

.122835 

.2 

0 

0 

0 

.013333 

.013333 

.040000 

.026667 

.066667 

.034477 

.107750 

a 

0 

0 

-.000417 

.002500 

.006251 

.026667 

.020000 

.053333 

.028476 

.092665 

.4 

0 

0 

0 

0 

0 

.013333 

.013333 

.040000 

.022476 

.077580 

a 

0 

0 

0 

0 

-.000417 

.002500 

.006251 

.026667 

.016477 

.062495 

.6 

0 

0 

0 

0 

0 

0 

0 

.013333 

.010476 

.047410 

.7 

0 

0 

0 

0 

0 

0 

-.000417 

.002500 

.004060 

.032325 

.8 

0 

0 

0 

0 

0 

0 

0 

0 

-.001523 

.017240 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

-.001077 

.  .004190 

TABLE  III.-  VALUES  OF  THE  INTEGRATING-  MATRIX 


NAG  A  TN  No.  1876 


(a)  Six— Point  Solution 


.08333  .13333 

-.01667 

0 

0 

06667  .26667 

.06667 

0 

0 

06667  .26667 

.15000 

.13333 

-.01667 

06667  .26667 

.13333 

.26667 

.06667 

06667  .26667 

.13333 

.26667 

.10833 

(b)  Ten— Point  Solution 


t)/sa  0 


.2 

.3 

0 

0 

-.00833 

0 

.03333 

0 

.07500 

.06667 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.06667 

.13333 

.5 

.6 

.7 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.06667 

-.00833 

0 

0 

0 

.13333 

.03333 

0 

0 

0 

.13333 

.07500 

.06667 

-.00833 

0 

.13333 

.06667 

.13333 

.03333 

0 

.13333 

.06667 

.13333 

.07500 

.06667 

\NACA^ 
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TABLE  IV.-  VALUES  OF  THE  LOAD-CONVERSION  MATRIX  Qcg] 


(a)  Six— Point  Solution 


i/sa 

0 

'  .2 

.4 

.6 

.8 

.9 

0.2 

0.08333 

0.13333 

-0.01667 

0 

0 

0 

.4 

-.01667 

.13333 

.08333 

0 

0 

0 

.6 

0 

0 

.08333 

.13333 

-.01667 

0 

.8 

0 

0 

-.01667 

.13333 

•08333 

0 

•9 

0 

0 

0 

0 

.04553 

.05752 

1.0 

0 

0 

0 

0 

-.01886 

.09333 

(t )  Ten— Point  Solution 


VSA 

0 

.1 

.2 

•3 

.4 

.5 

.6 

.7 

.8 

.9 

0.1 

0.04166 

0.06667 

-0.00833 

0 

0 

0 

0 

0 

0 

0 

.2 

-.00833 

.04167 

0 

0 

0 

0 

0 

0 

0 

.3 

0 

.04166 

►06667 

-.00833 

0 

0 

0 

0 

0 

0 

-.00833 

c 

06667 

.04167 

0 

0 

0 

0 

0 

.5 

0 

0 

0 

.04166 

. 

06667 

-.00833 

0 

0 

0 

.6 

0 

0 

0 

-.00833 

• 

06667 

.04167 

0 

0 

0 

.7 

0 

0 

0 

0 

0 

.04166 

.06667 

-.00833 

0 

.8 

0 

0 

0 

0 

0 

0 

-.00833 

.06667 

.04167 

0 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

.04166 

0 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

-.00833 

.09333 
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TABLE  V.—  VALUES  OF  THE  LOAD-CONVERSION  MATRIX 


(a)  Six— Point  Solution 


0.2  0.01667  0.16667  0.01667  0  0  0 

.4  -.00833  .05000  \  11667  .05000  00833  0 

.6  0  0  .01667  .16667  .01667  0 

.8  0  0  -.00833  .05000  .08946 

.90  0  0  0  .00631 

1.0  0  0  0  0 _ — .  01077 - 


(t>)  Ten— Point  Solution 


0 

cl 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

0.1 

0.00833! 

0.08333 

b.  00831 

0 

0  1 

■  I 

0 

0 

0 

0 

0 

.2 

-.00417 

002500 

.05834 

.02500 

c- 

H 

0 

0 

• 

1 

0 

■  1 

0 

0 

0 

0 

H 

0 

0 

.00833 

.08333 

.00833 

0  1 

; 

0 

0 

0 

0 

D 

0 

0 

-.00417 

.02500 

.05834 

.02500 

-.00417 

0 

0 

0 

•5 

0 

0 

0 

0 

.00833 

.08333 

.00833 

0 

0 

0 

.6 

0 

0 

0 

0 

-.00417 

.02500 

.05834 

.02500 

-.00417 

0 

.7 

0 

0 

0 

0 

0 

0 

.00833 

.08333 

.00835 

0 

.8 

0 

0 

0 

0 

0 

0 

-.00417 

.02500 

.06029 

.02035 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

.08860 

1.0 

0 

0 

0 

0 

0 

0 

0 

0 

-.01077 

.04190 
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TABLE  VI.-  FORM  FOR  COMPUTATION  OF  AEROELASTIC  MATRIX  -  Concluded 


(b)  Computing  Instructions  —  Concluded 


Supersonic  Case 


h2 


NACA  TN  No.  1876 

ft 

TABLE  VII.-  FORM  FOR  SOLUTION  OF  AEROELASTIC  EQUATION 


(L)  Aerodynamic  Loading 


(a)  Divergence 


Auxiliary  matrices 


IT  AC  A  TN  Wo.  1876 


TABLE  IT.-  COMPUTATION  OF  AEROELASTIC  MATRIX  OF  EXAMPLE  WING 
(SUBSONIC  CASE) 


© 

1  f(GJ)r 

w 

© 

eirvW 

>“ 

_sub 

fl/8* 

0 

.2 

.4 

.6 

.8 

.9 

l/sA 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

0.06667 

0.21494 

0.08560 

0.13120 

0.03435 

0.04707 

.2 

.08333 

.15066 

-.03501 

0 

0 

0 

.2 

-.01667 

.10746 

.09630 

.13120 

.03435 

.04707 

.06667 

.30134 

.14001 

0 

0 

0 

.4 

0 

0 

.04280 

.13120 

.03435 

.04707 

.6 

.06667 

.30134 

.31500 

.58665 

-.16070 

0 

.6 

0 

0 

-.01070 

.06560 

.04048 

.04707 

.8 

.06667 

.30134 

.27999 

1.17335 

.64270 

0 

.8 

0 

0 

0 

0 

.00981 

.04707 

•9 

.06667 

.30134 

.27999 

1.17335 

1.04430 

1.04072 

.9 

0 

0 

0 

0 

-.00694 

.02912 

© 

L1 

1  Red 

'3j  |_EI 

© 

©] 

[d 

)] 

i]/sA 

- : - 

0 

.2 

.4 

.6 

.8 

.9 

n/sA 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.2 

.08333 

.14933 

-.03201 

0 

0 

0 

.2 

-.00613 

— 

.01777 

.01728 

.01922 

.00503 

.00689 

.4 

.06667 

.29867 

.12801 

0 

0 

0 

.4 

-.01173 

_ 

.03724 

.04875 

.08089 

.02118 

.02902 

.6 

.06667 

.29867 

.28800 

.49465 

-.12452 

0 

.6 

-.01173 

.03724 

.05082 

.17174 

.06273 

.07104 

.8 

.06667 

.29867 

.25599 

.98935 

.49802 

0 

.8 

-.01173 

.03724 

.03938 

.22118 

.10693 

.16293 

•9 

.06667 

.29867 

.25599 

.98935 

.80923 

.73070 

.9 

-.01173 

.03724 

-.03938 

.22118 

.10254 

.23190 

® 

|s^(tan2A) 

© 

MW. 

l/flA 

0 

.2 

.4 

.6 

.8 

.9 

4/sA 

0 

.2 

.4 

.6 

.8 

.9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.04826 

0.04325 

0.11456 

0.04549 

0.07917 

.2 

.04585 

.08216 

-.01761 

0 

0 

0 

.2 

-.00167 

.00905 

.02028 

.07638 

.03389 

.06185 

.4 

.03668 

.16433 

.07043 

0 

0 

0 

.4 

0 . 

0 

0 

.03818 

.02227 

.04453' 

.6 

.03668 

.16433 

.15846 

.27216 

-.06851 

0 

.6 

0 

0 

-.00135 

.00716 

.00963 

.02721 

.8 

.03668 

.16433 

.14085 

.54434 

.27401 

0 

.8 

0 

0 

0 

0 

-.00095 

.00990 

•9 

.03668 

.16433 

.14085 

.54434 

.44524 

.40203 

.9 

0 

0 

0 

0 

-.00067 

.00241 

@ 

[©][ 

©j 

*l/BA 

0 

.2 

.4 

.6 

.8 

.9 

© 

tan  A)  L*°J 

0 

0 

0 

0 

0 

0 

0 

n/sA 

0 

.2 

.4 

.6 

.8 

.9 

.2 

-.00199 

.01106 

.02444 

.08202 

.03500 

.06286 

0 

0 

0 

0 

0 

0 

0 

.4 

-.00399 

.01543 

.04289 

.20681 

.09790 

.18314 

.2 

-.1629 

0 

0 

0 

0 

0 

.6 

-.00399 

.01543 

.03755 

.28397 

.16542 

.33786 

.4 

-.1629 

0 

0 

0 

0 

0 

.8 

-.00399 

.01543 

.03221 

.30252 

.19308 

.48336 

.6 

-.1629 

0 

0 

0 

0 

0 

.9 

-.00399 

.01543 

.03221 

.30252 

.18680 

.52208 

.8 

-.1629 

0 

0 

0 

0 

0 

[© 

DO  -  f© 

]-[©] 

•9 

-.1629 

0 

0 

0 

0 

0 

n/sA 

0 

.2 

.4 

.6 

.8 

.9 

&  (<W  -  1»aMtan  A)  "  -°'lffl9  ■ 

0 

0 

0 

0 

0 

0 

0 

© 

©>[©>[©] 

.2 

-.00414 

.00671 

-.00716 

-.06280 

-.02997 

-.05597 

*l/»A 

0 

.2 

.4 

.6 

.8 

.9 

.4 

-.00774 

.02181 

.00586 

-♦12592 

-.07672 

-.15412 

0 

0 

0 

0 

0 

0 

0 

.6 

-.00774 

.02181 

.01327 

-.11223 

-.10269 

-.26682 

.2 

-.03372 

.23282 

-.05262 

0 

0 

0 

-.00774 

.02181 

.00717 

-.08134 

-.08615 

-.32043 

•4 

-.05955 

.46567 

.21044 

0 

0 

0 

•9 

-.00774 

.02181 

.00717 

-.08134 

-.08426 

-.29018 

.6 

-.05955 

.46567 

.47346 

.85881 

-.22921 

0 

.8 

-.05955 

.46567 

.42084 

1.71769 

.91671 

0 

.9 

-.05955 

.46567 

‘.42084 

1.71769 

1.48954 

1.44275 
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TABLE  I.-  SOLUTION  OF  AEROELASTIC  EQUATION  FOE  EXAMPLE  WING  (SUBSONIC  CASE) 


(b)  Aerodynamic  Loading 


(a)  Divergence 


4--°*8s  a“0-552 


[W- 

.[a]] 

Tl/sA 

0 

.2 

.4 

.6 

.8 

•9 

0 

l.oooo 

0 

0 

0 

0 

0 

m 

.0023 

.9963 

.0040 

.0347 

.0165 

.0309 

.4 

.0043 

.9968 

.0424 

.0851 

.6 

.0043 

-.0073 

1.0620 

.0567 

.1473 

.8 

.0043 

-.0120 

m 

.0449 

1.0476 

.1769 

.9 

.0043 

-.0120 

.0449 

.0465 

1.1602 

0 

D 

.4 

.6 

.8 

•9 

• 

0 

0 

0 

0 

0 

0 

n 

-.0041 

.0067 

J 

-.0300 

-.0560 

D 

-.0077 

.0218 

.0059 

-.0767 

-.1541 

EJ 

-.0077 

.0218 

.6133 

-.1027 

-.2 668 

.8 

-.0077 

.0218 

.0072 

-.0862 

-.3204 

m 

-.0077 

.0218 

.0072 

MB 

-.2902 

R 

N» 

N- 

1 

0.1033 

1 

.2826 

1 

.4620 

1 

.6413 

1 

.8207 

1 

.9103 

Auxiliary  matrices 


H 

(1) 

(2) 

(3) 

(4) 

(5) 

0 

0 

0 

0 

0 

0 

a 

.3000 

.3115 

.3449 

.3480 

.5000 

.7311 

.7846 

.7876 

HI 

m 

.7000 

1.0286 

1.0526 

1.0532 

■H 

a 

.9000 

1.0775 

1.0714 

1.0713 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

a 

1.0000 

0 

0 

0 

0 

0 

B 

.0023 

.9963 

.0040 

.0348 

.0166 

.0310 

B 

.0043 

.9968 

.0702 

.0427 

.0857 

B 

.0043 

5§gf 

1.0629 

.0538 

.1395 

B 

.0043 

-.0120 

.0039 

.0456 

1.0455 

.1638 

B 

.0043 

-.0120 

.0039 

.0456 

.0444 

1.1473 

1.0000 

0.1033 

1.0014 

.2834 

1.0110 

.4665 

.9551 

.6094 

.9261 

.7630 

.8081 

.7439 

["■s  -1 

Ra}a 

Pe  J 

K  Ip 

J“e  -1 

Kt 

1.0000 

0.1033 

.9320 

.2320 

.8518 

U> 

£ 

ut 

.7996 

.4711 

.7937 

.6412 

.8081 

.7439 

N 

|i! 

.  0 
c/°r] 

mmi 

0.1909 

0.0581 
. - 

0.0866 

I 

A1W 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

0 
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letricai  parameters  used  in  the  analysis. 


Rotations  due  to  root 
deformations,  fr  and  rr, radians 


HACA  TN  Ho.  1876 


I 


104  in. 
w=  15  in. 

(GJ\  =  5.10  10®  lb-in.1 
3.47'*  10  6  lb-in.! 

Tr  =  43,470  in- lb 
Mr  =  260,000  in.- lb 


Rotation 

Av.  value 

Q  val 

?rT 

-0.00 02 

0.16 

Pru 

.0025 

.33 

— 

—  -  rrT 

.0010 

1.45 

— 

-  r rM 

-  .0010 

-  .24 

Distance  along  span,  17  ,  inches 

Figure  2..-  Rotations  of  a  45°  swept  box  beam  due  to  root 
deflect  ion  5  (data  from  reference  5). 
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0  .2  .4  .6  .8  1.0 

Lateral  ordinate,  y/sw 
(a)  Assumed  rigid -wing  loading. 
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Lateral  ordinate,  y/sw 


(b)  Flexible -wing  loading. 

Figure  5.-  Loqd  distribution  of  example  wing- 


